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I.   INTRODUCTION 

A.  NEED  FOR  STUDY 

Electromagnetic  pulse  has  become  a  possible  threat  to 
nearly  all  sophisticated  military  systems  [Ref.  1].   Under 
the  proper  circumstances  a  significant  portion  of  the  energy 
released  during  a  nuclear  detonation  will  be  appear  as  an 
Electromagnetic  Pulse  (EMP)  having  the  same  frequencies  as 
those  employed  by  most  of  our  commercial  radio  and  military 
system  equipments. 

Any  metallic  object,  such  as  an  aircraft  structure,  ex- 
posed to  electromagnetic  fields  can  be  a  collector  of  elec- 
tromagnetic energy,  that  is,  act  like  an  antenna,  even 
though  it  was  never  intended  to  be  that. 

The  crossed-dipole  receiving  antenna  has  been  used  as 
a  representative  model  to  approximate  electromagnetic  pulse 
effects  on  aircraft  [Ref.  2].   The  crossed-dipole,  in  this 
case  the  aircraft,  can  be  modeled  in  its  simplest  way  as  a 
two-dimensional  structure,  that  is,  a  cross  with  thickness 
equal  to  zero. 

B.  BACKGROUND  INFORMATION 

Electromagnetic  boundary-value  problems  can  be  solved 
exactly  via  classical  separation  of  variable  analysis  only 
for  a  few  geometries,  such  as  the  sphere,  spheroid,  circular 
cylinder,  elliptical  cylinder,  strip  and  wedge.   Such  solu- 
tions can  be  expressed  in  terms  of  a  summation  of  a  set  of 
eigenfunctions  which  can  be  evaluated  with  a  high  speed 
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computer.   However,  for  problems  involving  complicated  geo- 
metries the  exact  solution  is  not  available.   Thus ,  the  ap- 
proximate and  numerical  methods  are  of  great  significance. 
The  research  discussed  herein  is  directed  toward  the  numeri- 
cal solution  of  the  problem  of  scattering  from  conducting 
bodies  of  arbitrary  shape. 

Two  methods  [Ref.  3]  are  available  for  electromagnetic 
modeling  of  continuous  conducting  surfaces  with  arbitrary 
shape:   the  wire  grid  model  and  the  surface  current  model 
using  rectangular-pulse  bases.   Both  methods  have  similar 
limitations  with  the  maximum  cell  width  restricted  to  ap- 
proximately A/10.   Unless  the  conducting  body  is  symmetric 
or  is  a  figure  of  revolution,  computer  storage  requirements 
have  limited  the  moment-method  application  to  bodies  with  . 
surface  area  not  exceeding  one  or  two  square  wavelengths. 

In  this  paper,  the  piecewise-sinusoidal  reaction  tech- 
nique applied  by  Richmond  [Ref.  3]  to  scattering  by  conduct- 
ing bodies  of  arbitrary  shape  is  used. 

C.   STATEMENT  OF  THE  PROBLEM 

The  objective  of  this  research  is  to  find  the  current 
distribution  in  a  two-dimensional  cross.   In  order  to  get 
an  understanding  of  the  "junction  problem"  that  will  be 
present  in  the  cross  structure,  the  induced  current  density 
on  a  perfectly-conducting  rectangular  plate  will  be  first 
solved  for,  this  result  having  been  presented  by  Richmond 
and  Wang  [Ref.  3],   Likewise  the  square  plate  solved  by 
Rahmat-Samii  and  Mitra  [Ref.  4]  will  be  studied.   From  this 
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the  complexity  of  the  structure  will  be  increased,  solving 
for  current  distribution  on  structures  like  an  L,  T  and 
finally  the  cross. 

D.   SCOPE  AND  LIMITATIONS  OF  THE  STUDY 
The  general  solutions  steps  are: 

1)  Divide  the  continuous  conducting  surface  into  cells 
with  the  surface  current  distribution  being  expanded  in 
overlapping  sinusoidal  bases  functions . 

2)  Via  an  application  of  Galerkin's  method,  the  inte- 
gral equation  formulated  with  the  zero-reaction  concept  is 
reduced  to  a  matrix  equation. 

3)  The  current  distribution  over  the  conducting  surface 
is  determined  via  matrix  inversion. 

The  manner  in  Wiiicii  t.ue  structure  is  subsectioneu.  is 
very  important  in  any  numerical  analysis  approach.   Likewise 
the  basic  current  distribution  chosen  for  the  subsections 
drastically  affects  the  amount  and  difficulty  of  calculations 
These  methods  will  yield  approximate  results,  the  quality  of 
the  approximation  depending  on  the  amount  of  effort  (and 
computer  time)  one  is  willing  to  devote  to  the  problem. 

As  will  be  seen  symmetry  may  simplify  considerably  the 
number  of  calculations  involved.   In  some  cases  one  may  want 
to  take  advantage  of  symmetry;  in  other  cases  it  may  be 
simpler  to  ignore  it.   It  appears  that  one  must  trade  some 
formulation  time  for  computer  time.   Is  the  tradeoff  worth- 
while?  The  answer  to  this  question  will  depend  on  the  extent 
of  the  computations.   In  reducing  the  number  of  unknowns,  one 
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increases  the  analysis  time  but  reduces  computer  programming 
time  and  computer  running  time. 
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II.   METHOD  OF  ANALYSIS 

A.  METHOD  OF  MOMENTS 

The  method  of  moments  has  been  widely  explained  in  the 
literature.   The  interested  reader  may  refer  on  this  subject 
to  References  5,  6  and  7. 

B.  PIECEWISE  SINUSOIDAL  BASIS  SETS 

The  physical  problem  may  be  represented  by  many  different 
operator  equations,  and  a  suitable  one  must  be  chosen.   There 
are  also  an  infinite  number  of  sets  of  expansion  functions 
and  testing  functions  that  may  be  chosen. 

Intuitively,  the  more  nearly  the  expansion  set  resembles 
the  actual  form  of  the  current  distribution,  the  fewer  the 
terms  required  to  accurately  represent  the  current.   In  ad- 
dition, certain  types  of  basis  functions  may  reduce  many 
time-consuming  numerical  operations  incorporated  in  approxi- 
mations.  The  piecewise-sinusoidal  basis  set  exhibits  these 
qualities.   The  objective  is  to  represent  the  current  dis- 
tribution, over  individual  subsections  of  a  domain  with  ex- 
pansion functions  that  possess  greater  curve-fitting  potential 
than  do  flat  pulses  and  yet  lead  to  well-conditioned  matrices 
in  approximate  solution  procedures.   The  application  of  piece- 
wise-sinusoids  in  Galerkin's  method  and  the  method  of  moments 
to  electromagnetic  field  problems  is  known  as  the  Piecewise- 
Sinusoidal  Reaction  Matching  Technique. 
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C.   PIECEWISE-SINUSOIDAL  REACTION  MATCHING  TECHNIQUE.  SURFACE  MODEL 

The  reaction  concept  and  its  application  has  been  dis- 
cussed by  Rumsey  [Ref.  8],  Cohen  [Ref.  9],  Harrington  [Ref. 
10]  and  Richmond  [Ref.  11]. 

The  following  material  has  been  previously  presented 
by  Richmond  and  Wang  [Ref.  3];  it  is  presented  here  for 
completeness . 

Consider  the  exterior  scattering  problem  illustrated 
in  Figure  1.   In  the  presence  of  a  dielectric  or  conducting 
body,  the  impressed  electric  and  magnetic  currents  (J. ,M . ) 
generate  the  electric  and  magnetic  field  intensities  (E  ,H) . 
For  simplicity  let  the  exterior  medium  by  free  space. 

From  the  surface-equivalence  theorem  of  Schelkunoff , 
the  interior  field  will  vanish  (without  disturbing  the  ex- 
terior field)  if  the  following  surface-current  densities 
are  introduced: 

J   =  n  x  H  (1) 

— s       — 

M   =  E  x  n  (2) 

— s    — 

on  the  closed  surface  s  of  the  scatterer.   (The  unit  vector 

n  is  directed  outward  on  s . )   In  this  situation,  illustrated 

in  Figure  2,  it  is  possible  to  replace  the  scatterer  with 

free  space  without  disturbing  the  field  anywhere. 

By  definition,  the  incident  field  (E. ,H. )  is  generated 

by  ( J. ,M . )  in  free  space,  and  the  scattered  field  is 

E   =  E  -  E.  (3) 

— s    —    —1 

H   =  H  -  H. .  (4) 

— s    —   — 1 
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Figure  1.   The  Source  (J . ,M. )  Generates  the  Field  (E,H)  with 
Scatterer . 


(0,0) 

SCATTERER  OR 
FREE  SPACE 


— s  — s 


(E,H) 


0 


Ji9M. 

^SOURCE 


Figure  2 .   The  Interior  Field  Vanishes  When  the  Currents 

(J  ,M  )  are  Introduced  on  the  Surface  of  the 
— s  - -s 
Scatterer. 
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When  the  surface  current  (J  ,M  )  radiates  in  free  space ,  it 

— s  — s 

generates  the  field  (E  ,H  )  in  the  exterior  and  C-E.,-H.) 
°  — s  — s  —1   —1 

in  the  interior  region.   This  result,  illustrated  in  Figure 
3,  is  deduced  from  Figure  2  and  the  superposition  theorem. 

With  the  scatterer  replaced  by  free  space,  it  can  be 
noted  in  Figure  2  that  the  interior  region  has  a  null  field. 
As  shown  in  Figure  M- ,  by  placing  an  electric  test  source  J 
in  this  region  it  is  found  from  the  reciprocity  theorem 
that 

#s(Is'Et-Ms-Ht)ds  +  ///(Ji*Et-Mi-Ht)dv  =  0        (5) 

where  (E ,H )  is  The  free  space  field  of  the  test  source. 
In  words,  Equation  (5)  states  that  the  interior  test  source 
has  zero  reaction  with  the  other  sources.   This  "zero  reac- 
tion theorem"  was  developed  by  Rumsey  [Ref.  8]. 

Equation  (5)  is  the  integral  equation  for  the  scattering 
problem,  and  the  objective  is  to  use  this  equation  to  determine 
the  surface  current  distribution  J   and  M  .   To  accomplish 
this,  these  functions  are  expanded  in  finite  series  yielding 
a  finite  number  M  of  unknown  expansion  constants.   Next  N 
simultaneous  linear  equations  are  obtained  to  permit  a  solu- 
tion for  these  constants.   One  such  equation  is  obtained  from 
Equation  (5)  each  time  a  new  test  source  is  set  up. 

If  the  scatterer  is  a  perfect  conductor,  the  magnetic 
current  M   vanishes.   Assuming  a  finite  conductivity  and 
using  the  impedance  boundary  condition 

M   =  Z   J   x  n  (6) 

— s    — s  — s 
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where  Z   denotes  the  surface  impedance  equation  (5)  yields 
s 

-#   J  «[E  -(nxH.  )Z  Ids  =  ///(J.  -E.  -M.  -H.  )dv.        (7) 

s   s   — t      — t   S  —1  — t  —1  — t 

For  two-dimensional  problems  involving  arbitrary  scat- 

terers ,  J   and  M   are  functions  only  of  the  position  on  the 
— s     — s 

surface  of  the  scatterer.   If  M.  vanishes  and  denoting 

(E  ,H  )  as  the  free-space  field  of  test-source  m,  Equation 
m'  m  r  '   H 

(7)  yields 

-//   J  -[E  -(nxH  )Z  ]ds  =  //  J.-Eds.  (8) 

s  — s   — m    — m   s  l  — m 

Now  representing  the  electric  current  distribution  as 

follows : 

N 

J   =   Z,  I   J  (9) 

— s    n=l   n  — n 

where  the  complex  constants  I   are.  samples  of  the  function 

J  .   The  vector  functions  J   are  known  as  basis  functions . 
— s  — n 

Basis  functions  J  ,  and  testing  (weighting)  sources  J    with 

— n '  &      &     &  _m 

unit  current  density  at  their  terminals  are  used. 

From  Equations  (8)  and  (9)  the  following  simultaneous 
linear  equation  is  obtained 

N 

£,  I   C    =  A  (10) 

n=l   n   mm    m 

with  m  =  1,2,3,. ..,N j  where 

C    =-//  J  .[E  -(nxH  )Z  ]ds 

mm      — n   — m    -m   s 
n 

=  -ff      J  -E   ds  (11) 

m  -m  — n 

A   =  //.  J.-E   ds  =  //  J  -E.  ds  (12) 

m     i  — i  — m        m  — m  — l 
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The  integrations  in  Equations  (11)  and  (12)  extend  over 

the  region  where  the  integrand  is  non-zero.   Region  "m" 

covers  the  interior  test  source,  J  ,  and  region  "n"  is  that 

— m        ° 

portion  of  the  surface  s  covered  by  the  expansion  function 

J  .   The  first  and  the  second  integrals  in  Equation  (11)  are 
— n 

related  by  the  reciprocity  theorem,  where  E   is  the  free- 
space  field  generated  by  J   and  the  associated  magnetic 

current  M  . 
— n 

D.   UNIDIRECTIONAL  CURRENT  MODEL 

A  two-dimensional  surface  dipole  located  on  the  y-z  plane 
is  illustrated  in  Figure  5(a).   This  source  is  an  electric 
surface-current  density  with  height  2D  and  width  w.   The 
surface  current  density  is  given  by 

T       sinK(D-z)   ~       ,.    n  ^    ^  n        /no> 

J,    x  =  s — tttt —   z       for  0<z<D        (13) 

— (z)      smKD 

T     _  sinK(D+z)   A       ,.    _  .    .  n        *,..» 

J,  x  =  s — ^ z       for-D  <  z  <  0        (14) 

— (z)      smKD 

where         K  =  2tt/A. 

As  illustrated  in  Figure  5(b),  the  current  density  is  as- 
sumed to  be  piecewise-sinusoidal  in  the  z-direction,  with 
period  equal  to  A,  vanishing  at  the  edges  z  =  D  and  z  =  -D, 
and  equal  to  unity  at  the  center.   From  Figure  5(c)  it  can 
be  seen  that  the  current  density  is  assumed  to  be  uniformly 
distributed  in  the  tranverse  direction. 

The  dipole  current  distribution,  Equations  (13)  and 

(14),  will  be  used  as  the  basis  functions  (J   in  Equation 

—n.  ^ 

(9))  for  expanding  the  unknown  current  distribution  induced 
on  a  conducting  surface.   Also,  surface  dipoles  will  be 
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Figure  5.   An  Electric  Surface  Dipole  and  Its  Current 
Density  Distribution. 


23 


employed  as  test  sources  with  the  reaction  concept  to  solve 
the  integral  equation. 

1 .   Calculation  of  the  Coupling  Term 

The  complex  number  C    in  Equation  (11)  represents 

c  mm      ^  c 

the  mutual  coupling  between  two  sources.   Let  <m,n>  denote 
the  mutual  coupling  between  sources,  then  Equation  (11)  can 
be  written  as 

C    =  -<m,n>  =  -ffJ    «E„  ds  (15) 

mn  m  — m  n 

Relating  the  mutual  coupling  between  two  sources  to  the  cir- 
cuit parameters  [Ref.  3] 

C    =  V    I  (16) 

mn     mn   on 

where  V    is  the  open  circuit  terminal  voltage  induced  at 
mn         r  & 

source  m  by  source  n,  and  I    is  the  short  circuit  terminal 
J  '       on 

current  induced  at  source  n. 

The  surface  dipole  self  impedance  has  been  defined  [Ref 
3]  with  the  induced-emf  formulation: 

V       C 

Z    =  J5E  =    ™    .  (17) 

mm    I      II 

om    om  om 

From  Equation  (15),  Equation  (17)  yields 

Z    =  ^^ //  J  -E   ds.  (18) 

mm    II     m  — m  — m 
om  om 

where  J   is  the  surface-current  density  of  source  m  and  E_ 
— m  J  —m. 

is  its  free  space  electric  field. 

The  mutual  impedance  between  two  surface  dipoles  is 
defined  [Ref.  3]  by 
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Zmn,  =  7 T ^  J   ,En  dS«  <19> 

mm    11     m  — m  — n 
om  on 

The  evaluation  of  the  integrals  in  Equations  (18) 
and  (19)  requires  use  of  numerical  integration  techniques. 
K.  Khemayodhin  has  developed  a  computer  program  to  evaluate 
these  integrals.   Details  and  results  are  shown  in  his 
Engineer's  Degree  Thesis  to  be  presented  to  the  Electrical 
Engineering  Department,  Naval  Postgraduate  School,  in 
December  19  74.   His  computer  program  has  been  used  in  this 
thesis  to  calculate  the  self  and  mutual  coupling  between 
surface  dipoles  . 

Table  I  lists  the  self  impedance  of  a  center-fed 
two-dimensional  dipole  as  a  function  of  size.   Table  II 
lists  the  mutual  impedance  between  the  two-dimensional  di- 
poles shown  in  Figure  7  as  a  function  of  their  relative 
position . 

From  Equations  (17)  and  (19)  the  mutual  coupling 
is  derived  as 


C    =1    I    Z   .  (20) 

mn     om  on   mn 

For  the  dipole  shown  in  Figure  5,1    is  the  total 

current  at  z  =  0 .   Since  J  =  1  at  the  center  of  the  dipole 

w 

I    =  /  J  dy  =  w.  (21) 

om   -/  —  J 
o 

For  identical  dipoles  and  from  Equation  (21),  Equa- 
tion (20)  yields 

C    =  w2  x  Z   .  (22) 

mn         mn 
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Table  I.   Self  Impedance  of  Center-Fed  Two-Dimensional  Dipole 
Shown  in  Figure  6 . 


w/X 

R(l,l) 

X(l,l) 

0.2 

7.91 

-  68.06 

0.3 

17.91 

-  33.08 

0.4 

32.47 

-  14.54 

0.5 

53.11 

-   4.67 

He 


w 


» 


Figure  6.   Two-Dimensional  Dipole 
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Table  II.   Mutual  Impedance  of  Center-Fed  Two-Dimensional 
Dipoles  Shown  in  Figure  7. 


y/X 


z/X 


0.0 


0;25 


0.50 


0.75 


0.0 


67. 3+J13.9 


38.9-J22.? 


-8.5-J26.6 


19.2+iU.? 


0.25 


53.3+J59.7 


29.4-J8. 6 


-9.U-J19.7 


-15.6+n6.4 


0.50 


23.5+J6. 3 


10.0-J9. 2 


-9 . 5-i  6  .  6 


-7.2  +  3B.9. 


0.75 


1.1-J7.9 


2.9-J5.9  [-5.9  +  U.7 0  .  9  t.  jJ_Q 


-»  y 
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H"  w  ~H 
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Figure  7.   Coupled  Two-Dimensional  Dipoles  with  D  =  0.25A 
and  w  =  0.25A. 
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2 .   Calculation  of  the  Excitation  Column 

In  the  matrix  equation  C    I   =  A  ,  the  complex 

^        mn  n    m'        ^ 

quantities  A   form  the  excitation  column.   Physically,  A 
is  the  reaction  between  the  impressed  source  and  the  surface 
dipole  m.   If  the  impressed  source  is  an  electric  line 
source  with  current  density  J.,  Equation  (12)  reduces  to 

Am  =  ff      J    -E.  ds.  (23) 

m     m  — m  — i 

Assuming  a  line  source  located  at  a  great  distance  from  the 

surface,  the  incident  field  may  be  regarded  as  a  plane  wave 

with 

E.  =  E   ejKx.  (24) 

—l    o 

For  E   =  1  and  at  x  =  0 ,  E .  =1.   For  the  two-dimensional 
o  '  —1 

dipole  shown  in  Figure  5(a),  Equation  (23)  reduces  to 


-f  f 


D 


A   =  J    /   J   dz  dy.  (25) 

m   J        J        — m 

X)     o 

Using  the  surface  current  density  assumed  in  Equation  (13) 
and  (14),  Equation  (25)  yields 

A   =  v    2.W  vn  (1  -  cosKD).  (26) 

m    K  sin  KD 

E.   GENERAL  CURRENT  MODEL 

A  two-dimensional  surface  quadrupole  located  on  the  y-z 
plane  is  illustrated  in  Figure  8(a).   This  source  is  an  elec- 
tric surface-current  density  with  height  2D  and  width  2  w. 
The  surface  current  density  is  given  by 

T      _  Sin  K(D-Z)   A       r:  r>     ^  s-     r\  fnn\ 

J,    x  =  1 £r- —  z     for     0  <  z  <  D  (27) 

— (z)      sin  KD 
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Figure  8.   A  Two-Dimensional  Electric  Quadrupole  and  its 
Current  Density  Distribution. 
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T       sin  K(D+z)  A      ~  n  .    ^  n       /««% 

^(z)  =   sin  KB z      for  -D  <  z  <  0       (28) 


T     _  sin  K(w-y)  ~      ^     n  .  /*«% 

J,  ,=  — : i7 — ^—  y     for   0<y<w      (29) 

-(y)     sin  Kw     ;  J 


T  _  sin  K(w+y)  *      -  _  _       /on>> 

J,  v  =  1 t? — *—  y      for  -w  <  y  <  0       (30) 

— (y)     sin  K  w    ^  J 

where  K  =  2ir/X. 

As  illustrated  in  Figure  8(b)  and  (c),  the  current  den- 
sity is  assumed  to  be  piecewise-sinusoidal  in  the  y-direction 
and  z-direction,  with  period  equal  to  A,  vanishing  at  the 
edges  and  equal  to  unity  at  the  center.   Comparing  this  struc- 
ture with  the  one  used  in  the  Unidirectional  Current  Model 
it  can  be  appreciated  that  the  quadrupole  is  a  combination 
of  the  dipoles  used  in  the  former.   Its  current  distribution 
will  be  the  product  of  the  two  functions  shown  in  Figures 
8(b)  and  8(c). 

K.  Khemaydhin  has  also  developed  a  computer  program  to 
evaluate  the  current  distribution  on  conducting  bodies  of 
arbitrary  shape.   Details  are  contained  in  his  aforementioned 
thesis.   This  computer  program  has  been  used  in  this  thesis 
to  calculate  the  current  distribution  in  two-dimensional 
structures  when  current  in  two  directions  was  assumed. 
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III.   COMPUTER  PROGRAM.   DISCUSSION,  DESCRIPTION  AND  RESULTS 

This  section  is  concerned  with  the  discussion  and  de- 
scription of  the  computer  programs  used  to  find  the  current 
distribution  in  the  several  two-dimensional  structures  to 
be  treated.   There  are  a  number  of  refinements  which  can  be 
introduced  to  improve  the  accuracy  of  the  solution,  as  well 
as  to  increase  the  effectiveness  of  the  use  of  the  computer 
memory.   The  principles  involved  are  simple  in  concept;  in 
each  case  an  integral  equation  is  replaced  with  a  matrix 
equation,  which  is  inverted  to  solve  for  the  unknown  source 
function . 

A.   UNIDIRECTIONAL  CURRENT  MODEL 

In  the  solution  for  current  distribution  in  the  several 
two-dimensional  structures  in  which  only  one  direction  of 
current  was  allowed,  it  was  assumed  that  a  line  source  is 
located  at  a  great  distance  from  the  surface  so  that  the 
incident  field  may  be  regarded  as  a  plane  wave.   The  plane 
wave  is  assumed  to  be  incident  at  broadside  and  in  the 
z-direction . 

The  number  of  subsections  is,  of  course,  a  very  important 
quantity  in  any  Numerical  Analysis  Approach.   One  of  the 
easiest  methods  of  improving  accuracy  is  merely  to  increase 
the  number  of  subsections ,  but  increasing  the  number  of  sub- 
sections increases  the  number  of  unknowns ,  the  time  of  com- 
puter programming,  the  computer  running  time  and  the  computer 
memory  space.   Thus  it  appears  that  it  is  a  tradeoff  between 
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the  extent  of  computations  and  the  accuracy  of  the  solution 
desired. 

For  all  the  structures  solved,  the  surface  was  divided 
into  an  integer  number  of  equal  two-dimensional  dipoles. 
Although  this  method  of  segmentation  gives  an  equal  resolu- 
tion at  the  sharp  edges  (where  the  current  density  is  vary- 
ing very  rapidly)  as  well  as  in  the  center  of  the  surface 
(where  the  current  density  is  almost  constant),  it  allows 
the  use  of  symmetry  for  the  calculations  of  the  coupling 
terms.   The  use  of  symmetry,  as  explained  in  each  of  the 
particular  cases,  has  saved  considerable  computational  time. 

1.   The  Rectangular  Plate 

The  first  problem  solved  was  the  rectangular  plate, 
for  which  Richmond  and  Wans  CRef.  3D  have  previous  results. 
The  rectangular  plate,  as  shown  in  Figure  9,  was  divided 
into  ninety  overlapping  dipoles.   The  current  density  dis- 
tribution was  found  with  Computer  Program  1  listed  in  Ap- 
pendix A.   Its  description  is  as  follows:   the  program 
written  for  a  rectangular  shaped  structure  is  a  general  pro- 
gram applicable  for  different  methods  of  segmentation.   The 
only  restriction  is  that  the  segments  (dipoles)  must  be 
equal.   The  number  of  dipoles  is  limited  by  the  size  of  the 
matrix  which  the  inverting  matrix  subroutine  can  handle, 
usually  10  0.   Subroutine  CMTRIN,  modified  to  handle  matrices 
up  to  200  by  200,  was  used  in  this  program. 

The  way  to  describe  the  rectangular  plate  is  by  en- 
tering the  following  parameters:  (a)  LI  is  the  number  of 
dipoles  in  a  row  (30  in  Figure  9);  (b)  L3  is  the  number  of 
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rows  of  dipoles  ( 3  in  Figure  9);  (c)  L7  is  the  total  number 
of  dipoles  in  the  surface  (90  in  Figure  9);  (d)  D  is  the 
half-length  of  the  dipole  (see  Figure  5(a));  (e)  w  is  the 
width  of  the  dipole  (see  Figure  5(a)). 

DO  LOOPS  100,  200,  300,  and  400  generate  the  coor- 
dinates of  the  source  dipole  and  the  test  dipole  as  shown 
in  Figure  7.   Each  dipole  is  at  one  time  a  test  dipole, 
with  all  the  other  dipoles  becoming  source  dipoles.   In 
calculating  the  self  impedance,  the  Complex  Function  ZSDPAR 
provides  a  0.001X  separation  between  dipoles  and  calculates 
this  closely-coupled  case  of  mutual  impedance  as  the  self 
impedance . 

The  number  of  times  the  calculation  of  the  coupling 
between  dipoles  has  to  be  performed  depends  on  the  number 
of  dipoles  on  the  surface.   The  time  each  calculation  of 
the  coupling  term  takes  in  the  IBM  360/67  computer  has  been 
estimated  to  be  an  average  of  1.6  seconds.   To  perform  the 
8100  calculations  for  the  rectangular  plate  would  take  216 
minutes  of  computer  time.   However,  noticing  that  the  coup- 
ling between  dipoles  depends  only  on  their  relative  posi- 
tion, it  is  possible  to  apply  some  symmetry  to  the  problem. 
Calculating  only  the  first  row  of  the  coupling  matrix,  that 
is,  the  coupling  between  dipole  1  and  all  other  dipoles  in 
the  surface  will  yield  sufficient  terms  for  filling  the  rest 
of  the  matrix.   So,  only  90  calculations  of  the  coupling 
term  need  to  be  performed,  reducing  the  computer  time  dras- 
tically to  2.4  minutes.   DO  LOOPS  450  and  466  perform  the 
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symmetry  comparison  by  comparing  the  relative  position  of 
the  dipoles . 

The  calculation  of  the  right  hand  side  vector  of 
Equation  (10)  is  performed  with  the  closed  form  of  Equation 
(26). 

The  coefficient  matrix  is  then  inverted  using  sub- 
routine CMTRIN. 

DO  LOOPS  2400  and  2500  perform  the  multiplication  of 
the  inverted  matrix  [C(m,n)]    by  the  excitation  column 
CA(m)]. 

DO  LOOP  3000  performs  the  calculation  of  the  current 
density  at  each  dipole  and  normalizes  it  to  the  Magnetic 
Field  Intensity 

E 

n  =  g2-  .  (3D 

y 

From  Equation  (24) ,  E   =  1,  then 

H   =  -  .  (32) 

y   n 

The  absolute  value  of  the  current  density  is  plotted 
using  the  three-dimensional  plot  sub-program  PLT3D1  and  is 
shown  in  Figure  10 . 

2.   The  Square  Plate 

The  two-dimensional  square  plate  is  shown  in  Figure 
11.   The  structure  was  divided  in  90  overlapping  dipoles. 

The  current  density  distribution  was  found  using 
Computer  Program  1  following  exactly  the  same  procedure  as 
in  the  rectangular  plate.   The  only  difference  is  in  the 
size  of  the  dipoles.   For  the  square  plate  D  =  0.25A  and 
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Figure  10.   Rectangular  Plate  (1A  by  0.5A)  Magnitude  of 
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Figure  11.   Square  Plate  Segmented  in  90  Overlapping  Di- 
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w  =  1/30X,  instead  of  D  =  0.125A  and  w  =  1/30X  as  used  in 
the  rectangular  plate. 

The  current  density  distribution  is  shown  in  Figure 
12  and  was  plotted  using  the  three-dimensional  plot  sub- 
program PLT3D1. 

3 .   'L' -Shaped  Structure 

The  ' L1 -shaped  structure  is  shown  in  Figure  13.   The 
structure  was  divided  into  60  dipoles. 

Computer  Program  2  listed  in  Appendix  A  was  used  to 
solve  the  current  density  distribution  for  the  fL' -shaped 
structure .   This  is  a  general  program  applicable  for  dif- 
ferent ways  of  segmentation  of  an  ' L' -shaped  structure.   It 
has  the  same  restrictions  as  Computer  Program  1  (rectangular 
plate) ,  in  fact  both  programs  are  quite  similar.   Computer 
Program  2  assumes  first  that  the  surface  is  a  squared  sur- 
face and  following  the  same  procedure  described  in  Computer 
Program  1,  calculates  the  coupling  between  dipoles  of  the 
squared  surface.   The  program  has  been  designed  in  this  way 
in  order  to  be  able  to  apply  symmetry. 

DO  LOOPS  501  through  507  assign  the  corresponding 
coupling  from  the  squared  surface  to  the  ' L' -shaped  surface. 
Then,  the  same  procedure  as  in  Computer  Program  1  is  applied 
to  calculate  the  right  hand  side  vector,  to  invert  the  co- 
efficient matrix,  to  multiply  the  inverted  matrix  by  the 
right  hand  side  vector  and  finally  to  calculate  the  current 
density  at  each  dipole. 

The  magnitude  of  the  induced  current  density  is 

plotted  using  the  three-dimensional  plot  sub-program  PLT3D1 

and  is  shown  in  Figure  14. 
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Figure  12.   Square  Plate  (IX  by  IX)  Magnitude  of  Induced 
Current . 
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Figure  13.   ' L1 -Shaped  Structure  Segmented  in  60  Overlapping 
dipoles.   Plane  Wave  Incident  at  Broadside. 
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Figure  14.   '  L ' -Shaped  Structure  (IX  by  IX )  Magnitude  of 
Induced  Current. 
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4.  The  "I"  -Shaped  Structure 

The  'T* -shaped  structure  was  segmented  in  64  dipoles 
as  shown  in  Figure  15. 

Computer  Program  3  listed  in  Appendix  A  was  used  to 
solve  the  current  density  distribution  for  the  'T' -shaped 
structure.   The  procedure  is  exactly  the  same  that  the  one 
used  in  the  ' L' -shaped  structure.   The  differences  are  in 
the  number  of  dipoles  used  and  the  width  of  the  dipole 
(1/32  X    instead  of  1/30  X    used  in  the  'L'). 

The  current  density  distribution  is  shown  in  Figure 
16  and  was  plotted  using  the  three-dimensional  plot  sub- 
program PL3TD1. 

5 .  The  Two-Dimensional  Cross 

The  two-dimensional  cross  was  segmented  in  80  over- 
lapping dipoles  as  shown  in  Figure  17.   Computer  Program  4 
listed  in  Appendix  A  was  used  to  solve  for  its  current  den- 
sity distribution.   The  procedure  is  exactly  the  same  that 
the  one  used  for  the  'L' -shaped  structure.   The  differences 
are  in  the  number  of  dipoles  used  and  the  size  of  the  di- 
poles.  For  the  two-dimensional  cross  D  =  1/8  X  and  w  = 
1/16  A. 

The  current  density  distribution  is  shown  in  Figure 
18  and  was  plotted  using  the  three-dimensional  plot  sub- 
program PL3TD1. 

B.   GENERAL  CURRENT  MODEL 

For  the  solution  of  current  distribution  in  the  several 
two-dimensional  structures  in  which  two  directions  of  current 
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Figure  15.   "T'-Sh 


T  -Shaped  Structure  Segmented  in  64  Overlapping 
Dipoles.   Plane  Wave  Incident  at  Broadside. 
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Figure  16.   'T' -Shaped  Structure  (IX  by  IX)  Magnitude  of 
Induced  Current . 
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Figure  17.   Two -Dimensional  Cross  Segmented  in  80  Overlapping 
Dipoles.   Plane  Wave  Incident  at  Broadside. 
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Figure  18. 


Two-Dimensional  Cross  Structure  Magnitude 
of  Induced  Current. 
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was  allowed,  the  Computer  Program  written  by  K.  Khemayodhin 
was  used.   This  thesis  was  concerned  only  with  the  design 
of  the  subroutines  to  generate  the  geometry  of  the  structures. 
For  all  the  solutions  an  electric  field  of  120tt  volts/meter 
incident  at  broadside  was  assumed. 

Computer  Program  5  as  listed  in  Appendix  A  was  used  to 
generate  the  geometry  of  all  the  structures.   This  program 
generates  a  square  plate  of  dimensions  1  X  by  1  X.   The 
difference  in  the  particular  structures  that  has  been  solved 
is  the  current  path  assumed.   Thus,  this  program  can  be  also 
used  for  the  Unidirectional  Current  Model  by  only  specifying 
currents  in  one  direction.   The  reason  for  not  using  it  in 
the  Unidirectional  Current  Model  is  the  much  longer  computer 
running  time  it  takes ,  since  symmetry  considerations  are  not 
fully  applied. 

All  the  results  have  been  plotted  using  the  three-dimen- 
sional plot  subprogram  PL3TD1. 

1 .  The  Rectangular  Plate 

The  structures  of  dimensions  1  X  by  0 . 5  X  was  divided, 
as  shown  in  Figure  19,  in  32  surface  (monopoles).   Eight  cur- 
rent paths  were  assumed  in  the  z-direction  and  four  in  the 
y-direction.   Figure  20  shows  the  magnitude  of  the  induced 
current  in  the  direction  of  the  field.   Figure  21  shows  the 
magnitude  of  the  cross-polarized  induced  current. 

2.  The  Square  Plate 

The  structure  of  dimensions  1  X  by  1  X  was  divided 
in  64  monopoles,  as  shown  in  Figure  22.   Eight  current  paths 
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were  assumed  in  the  z-direction  and  eight  in  the  y-direction. 
Figure  2  3  shows  the  magnitude  of  the  induced  current  in  the 
direction  of  the  field.   Figure  2  4  shows  the  magnitude  of 
the  cross-polarized  induced  current. 
3.   The  'L' -Shaped  Structure 

The  ' L* -shaped  structure  was  divided  in  4  8  monopoles 
as  shown  in  Figure  25.   The  dimensions  of  the  horizontal  arm 
are  1  X  by  0.5  X,  and  the  vertical  arm  is  0.5  X  by  0.5  X. 
Eight  current  paths  were  assumed  in  the  z-direction  and 
eight  in  the  y-direction.   Figure  26  shows  the  magnitude  of 
the  induced  current  in  the  direction  of  the  field.   Figure 
27  shows  the  magnitude  of  the  cross-polarized  induced 
current . 

Li  ^"h  c      '  T  •      ov,-.r.-j     c-i — . ,  „-i-,  ,,~0 

H  .  ilic  x      -jilabcu     jLi'd^LuFc 

The  'T' -shaped  structure  was  divided  in  48  monopoles 
as  shown  in  Figure  28.   The  horizontal  arm  is  1  X   by  0 . 5  X 
and  the  vertical  arm  is  0.5  X    by  0.5  X.   Eight  current  paths 
were  assumed  in  the  z-direction  and  eight  in  the  y-direction. 
Figure  29  shows  the  magnitude  of  the  induced  current  in  the 
direction  of  the  field.   Figure  30  shows  the  magnitude  of 
the  cross-polarized  induced  current. 
5 .   The  Two-Dimensional  Cross 

The  two-dimensional  cross  was  divided  in  48  monopoles 
as  shown  in  Figure  31.   The  horizontal  arm  is  1  X  by  0.5  X 
and  the  two  vertical  arms  are  0.5  X  by  0.25  X.   Eight  current- 
path  were  assumed  in  the  z-direction  and  eight  in  the  y-di- 
rection.  Figure  32  shows  the  magnitude  of  the  current  induced 
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7.44 


Figure  23.   Square  Plate  (1  A  by  1  X).   Current  in  the 
Direction  of  the  Field. 
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Figure  24-.   Square  Plate  (1  X  by  1  X)  Cross-Polarized  Current 
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Figure  26.   'L'-Shaped  Structure  (IX   by  IX)  Current  in 
the  Direction  of  the  Field. 
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Figure  27.   'L' -Shaped  Structure  (1A   by  IX)  Cross-Polarized 
Current . 
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Figure  29. 


' T' -Shaped  Structure  (IX   by  IX)  Current  in  the 
Direction  of  the  Field. 
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Figure  30.   'T' -Shaped  Structure  (IX   by  IX)  Cross-Polarized 
Current . 
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Figure  32.   Two-Dimensional  Cross  (IX   by  IX)  Current  in 
the  Direction  of  the  Field. 
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in  the  direction  of  the  field.   Figure  33  shows  the  magnitude 
of  the  cross-polarized  induced  current. 
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Figure  33.  Two-Dimensional  Cross  (lx  by  lx)  Cross-Polarized 
Current. 
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IV.   NUMERICAL  RESULTS 

The  three-dimensional  plots  shown  in  the  preceding  section 
give  a  very  good  representation  of  the  current  density  dis- 
tribution over  the  two-dimensional  structure.   However,  due 
to  the  fact  that  some  of  the  points  plotted  come  from  inter- 
polation and  that  the  sub-program  PL3TD1  is  designed  to  plot 
only  rectangular  or  squared  surfaces,  some  accuracy  is  lost. 

In  the  tables  shown  in  this  section,  the  exact  values  at 
the  specific  points  are  given. 

For  the  Unidirectional  Current  Model  the  value  of  cur- 
rent is  at  the  center  of  the  dipole.   For  the  General  Cur- 
rent Method,  the  value  of  current  is  at  the  edges  of  the 
monopoles . 
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Table  III 


Rectangular  Plate  (Unidirectional  Current  Model) 
Magnitude  of  the  Induced  Current  at  the  Center 
of  the  Dipole.   See  Figure  9. 
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2.73                 1.57 

1.43 

2  .57 

1.43 

4.  38 
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4.38 

Table  IV.   Square  Plate  (Unidirectional  Current  Model). 

Magnitude  of  the  Induced  Current  at  the  Center 
of  the  Dipole.   See  Figure  11. 
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8.75     1   4.77 
1 1 

Table  V. 


' L' -Shaped  Structure  (Unidirectional  Circuit  Model) 
Magnitude  of  the  Induced  Current  at  the  Center  of 


the  Dipole.   See  Figure  13. 
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2.63        1. 39 
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3.71       29.92 

24.17 

5.61 

6.15 

' 

4.53 

3.60 

1 

5.56 

7.56 

33.51 

— — . — 

Table  VI. 


' T' -Shaped  Structure  (Unidirectional  Current 
Model).  Magnitude  of  the  Induced  Current  at 
the  Center  of  the  Dipole.   See  Figure  15. 
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9.30 

4.05 

4.78 
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3.16     U.20 
1 

3.16 

12.8     116.7     112.8 
1 

Table  VII. 


Two-Dimensional  Cross  (Unidirectional  Current 
Model).   Magnitude  of  Induced  Current  at  the 
Center  of  the  Dipole.   See  Figure  17. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  shown  for  the  rectangular  plate  (Unidirec- 
tional Current  Model)  has  been  compared  with  those  presented 
by  Richmond  and  Wang  [Ref.  3].   The  results  shown  for  the 
square  plate  (General  Current  Model)  has  been  compared  with 
those  presented  by  Rahmat-Samii  and  Mittra  [Ref.  4].   Both 
comparisons  are  seen  to  be  in  good  agreement. 

A  comparison  of  the  results  for  the  Unidirectional  Cur- 
rent Model  and  the  General  Current  Model  for  similar  geo- 
metries ,  suggest  that  a  solution  for  the  current  distribution 
in  two-dimensional  structures  is  not  complete  if  the  cross- 
polarized  current  is  neglected.   This  problem  is  more  re- 
markable where  the  junction  problem  is  present,  such  as  in 
the  ' L'  ,  TT'  and  cross  shaped  structures.   In  these  struc- 
tures the  magnitude  of  the  cross-polarized  current  reaches 
values  of  the  same  order  of  magnitude  as  that  for  the  cur- 
rent in  the  direction  of  the  field,  and  in  some  cases  as 
for  the  'T' -shaped  structure  and  the  two-dimensional  cross 
the  maximum  values  found  for  the  cross-polarized  current 
are  somewhat  larger  than  those  of  the  current  in  the  direc- 
tion of  the  field. 

The  computer  programs  used  in  the  General  Current  Model 
do  not  fully  apply  symmetry  in  the  calculation  of  the 
coupling  between  dipoles .   Since  this  calculation  is  the 
predominant  in  computer  time  consumption,  some  improvements 
should  be  investigated  in  order  to  be  able  to  solve  more 
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complicated  geometries  without  the  restriction  of  "computer 
time . " 

A  program  of  surface  current  measurement  for  flat-plate 
structures  with  corner  notches,  similar  to  the  fT!,  'L'  and 
cross  structures,  should  be  valuable  for  validating  these 
and  other  computer  codes  for  two-dimensional  surfaces. 
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APPENDIX  A:   FORTRAN  PROGRAM  LISTING 
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